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Using the density matrix formalism, we prove an existence theorem of the periodic steady-state 
for an arbitrary periodically-driven system. This state has the same period as the modulated 
external influence, and it is realized as an asymptotic solution (t—^--l-oo) due to relaxation processes. 
The presented derivation simultaneously contains a simple computational algorithm non-using both 
Floquet and Fourier theories, i.e. our method automatically guarantees a full account of all frequency 
components. The description is accompanied by the examples demonstrating a simplicity and high 
efficiency of our method. In particular, for three-level A-system we calculate the lineshape and 
field-induced shift of the dark resonance formed by the field with periodically modulated phase. 
For two-level atom we obtain the analytical expressions for signal of the direct frequency comb 
spectroscopy with rectangular light pulses. In this case it was shown the radical dependence of 
the spectroscopy lineshape on pulse area. Moreover, the existence of quasi-forbidden spectroscopic 
zones, in which the Ramsey fringes are signihcantly reduced, is found. The obtained results have a 
wide area of applications in the laser physics and spectroscopy, and they can stimulate the search 
of new excitation schemes for atomic clock. Also our results can be useful for many-body physics. 
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For the last several decades a rapid scientific-and- 
technological progress was substantially connected with 
an expansion of the lasers and laser technologies at dif¬ 
ferent areas of the science, engineering and industry. In 
this process the laser physics, nonlinear optics and laser 
spectroscopy take a special place. Many impressive suc¬ 
cesses in these directions are due to the theoretical sup¬ 
port, motivation and interpretation of experimental re¬ 
searches. In this context a paramount importance has 
the formulation of mathematical models (equations) and 
finding of their solutions, which adequately describe the 
physical picture of investigated problems. As an exam¬ 
ple, for atomic mediums the density matrix formalism 
is the most widespread approach describing the atom- 
field interaction and different relaxation processes (spon¬ 
taneous, collisional, etc.). Especial significance has so- 
called steady-state, which arises under the interaction of 
quantum system with stationary external fields. During 
long time this state plays a key role for theoretical de¬ 
scription of the basic problems in the laser physics and 
spectroscopy n 

However, in the last years the devices, in which dif¬ 
ferent parameters of electromagnetic fields are periodi¬ 
cally modulated, have gotten a great importance. First 
of all, it is so-called frequency comb-generators using 
the periodic pulse modulation of the laser field. Such 
sources of the pulse radiation are actively used now in the 
modern atomic clocks for frequency measurements an, 
and also they have promising per^ectives for the direct 
frequency comb spectroscopy a-Q. Also it is widely 
used now the phase (frequency) and/or amplitude pe- 
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riodic modulation of the laser field to form the reference 
resonances for small-size atomic clocks and magnetome¬ 
ters [IQ. Other promising and interesting direction in 
atomic clocks and magnetometers is connected wi th p eri- 
odically modulated polarization of the laser field . 

Of late years the understanding of periodically-driven 
systems is also one of the most active areas of research in 
many-body physics (e.g., see (2ll - l^ b In all above exam¬ 
ples the standard conception of steady-state is inapplica¬ 
ble, generally speaking, due to a regular time variation 
of external influence. 

In this paper, we generalize the steady-state concep¬ 
tion for an arbitrary quantum system under any periodic 
external influence. In this way we prove the following 
existence theorem; if the coefficients of density matrix 
dynamic equations have the period T, then the periodic 
solution with the same period T exists always. Due to 
the relaxation processes this solution is realized as an 
asymptotics (^—>-- 1 - 00 ) and therefore can be characterized 
as periodic steady-state. It is found the simple algorithm, 
which allows us directly to construct this solution inde¬ 
pendently of initial conditions and without the use of 
both Floquet and Fourier theories. Our approach con¬ 
siderably simplifies the analysis of the atom-field inter¬ 
action regardless of periodic modulation character: from 
smoothly harmonic type to the ultrashort pulses. As a 
bright example, for two-level atoms we analytically calcu¬ 
late the signal of the direct frequency comb spectroscopy 
formed by the periodic sequence of rectangular pulses. 
The obtained results have great perspectives in the laser 
physics and spectroscopy, and they can stimulate the 
search of new excitation schemes for atomic clock. 

The general mathematical formulation of the problem 
is following. For an arbitrary quantum system (e.g., 
atoms) let us consider the density matrix p(t) and its 
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normalization condition, which in the basis of states {I j)} Then, in accordance with Eq. Q, for other instant of 
can be presented as: time ^2 can write: 


Kt) =^\a)Pab{^){% Tr{p(t)} = = 1, (1) 

a,b j 


where p^bi^) s-re the matrix elements. The dynamic equa¬ 
tion for density matrix we reformulate to the vector form: 

dtp{t) = L{t) pit ); Tr{p(<)} = (n, pit)) = 1, (2) 

where the linear operator L(t) describes the interaction 
with external fields and different relaxation processes 
(spontaneous, collisional, etc.), the column-vector pit) 
is formed by the matrix elements Pabi^) using some deh- 
nite rule. As an example, for two-level system with basis 
states |1) and |2) there are four matrix elements p^bit) 
(a, 6=1,2), which can be ordered as following vector: 


P{t) 


( Piiit) \ 
Pl2(t) I 
P2l(t) I 

V 6*22(0 / 


(3) 


Besides we dehne the supplementary column-vector n 
permitting to express the value Tr{p(t)} as the dot prod¬ 
uct in Eq. where we use the standard dehnition of 
the dot product for arbitrary complex-valued vectors: 
(a:, y) = x'^y^. The vector n is formed in the follow¬ 

ing way: for the positions corresponding to the diagonal 
elements Pjjit) [as components of the vector pit)] the 
vector n has the value 1, and it equals to 0 for all other 
positions. In particular, for two-level system according 
to the definition ([3|) we obtain: 


n = 


(4) 


The normalization condition in @ implies the degen¬ 
eracy of equations [right-hand members, Lit)pit)] that 
corresponds to the following expression: 

(u,A(i)p)=0 (5) 


Ah) = Ah) ^ ( 8 ) 

where the two-time evolution operator A(t 2 ,ti) is deter¬ 
mined by the matrix Lit). Note that in the case of peri¬ 
odicity condition ([S]) we have the following relationship 

A(<2 + T, ti-|-T) = A(<2, ti) (9) 


for arbitrary t]^,t 2 - 

However, before the consideration of periodic case, let 
us prove several general statements. For this purpose, we 
multiply the Eq. ([5]) by the vector n as the dot product: 

{n,Ah)) = (n, Ait.2hi) Ah)) = [^Hhhi) fi, Ah)) > 

( 10 ) 

where AAt 2 ,ti) is the Hermitian conjugate operator 
to the A(t 2 ,ti). Because of (n, p(< 2 ))=(n, ^(ti)), from 
Eq. (fTUl) we obtain: 

in, Ah)) = {^Ahhi) A Ah)) ( 11 ) 

for arbitrary A^i)- This formula directly implies 

A\t2,thn = n, ( 12 ) 

i.e. the operator Al(t 2 ,ti) always has the eigenvector 
with the real eigenvalue 1. Therefore, the same eigen¬ 
value 1 also exists for the direct operator A(t 2 ,ti), i.e. 
always there is the eigenvector ^(^ 2 )^ 1 )) which satisfies 
the equation: 


MhhhAhhi) = Ahhi) ■ (i3) 

In contrast to Eq. m, the eigenvector r(< 2 ,ti) depends 
on t 2 , ti in the general case. 

Let us return to the case of periodicity dH). Consider 
the Ah at arbitrary instant of time t. In conformity with 
Eq. ([5]), the vector ph + T) is determined as 

At + T)=Ait + T,t)pit)- (14) 


for an arbitrary vector argument p. Note that the exis¬ 
tence of the time-independent vector ri is a base point for 
the next reasonings. 

Let us suppose an existence of the time period T in 
the operator Lit): 

Lit + T) = Lit). (6) 

In this case, as it will be shown below, the equation @ 
always has a periodic solution with the same period T: 

At+T) = m- ( 7 ) 

To begin the prove of this theorem, we assume that at 
some instant of time ti we have an arbitrary vector A^i)- 


Supposing the existence of the periodic solution pit + 
T)=pit), it follows from Eq. (IT4l) that this solution sat¬ 
isfies the equation: 

Ah = M't + T,t)Ahi (ii,p(i)) = i, (15) 

which always has nonzero solution due to the above 
proven statement of existence of eigenvector with eigen¬ 
value 1 for the operator A(t 2 ,ti) for arbitrary ti,t 2 [see 
Eq. (US])]- Using the Eqs. ® and dTTt . dTSl) . it can be 
easily shown that p(t -I- IT)—pit) (/=±1, ±2, ...). Tak¬ 
ing into account an arbitrariness of the time t in the 
Eqs. we can assert that the existence theorem 

of the periodic solution (0 is proven. 
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FIG. 1: The partition of the time interval [to, to + T] at 
subintervals and symbolic approximation of the dependence 
L{i) by the step function (red step line). 

Due to the relaxation processes there is an unique (in 
the majority of cases) periodic solution, which is real¬ 
ized during the time evolution as an asymptotic state 
(t^+oo) independently of initial conditions (similarly to 
the well-known case of stationary field (fT^ - see below). 
The solution © can be called as dynamic steady-state 
because of its dependence on t. 

However, taking into account a mathematical general¬ 
ity, one can theoretically assume the possibility of sev¬ 
eral solutions Pgit)} for Eq. (fT^ . when 

the eigenvalue 1 for the matrix A{t + T, t) is degenerated. 
In this case the general periodic solution has a form of 
superposition: 

Q Q 

p{t) ='^aqpq{t)-, {n,pq{t)) = l; ^ a, = 1,(16) 

q=l q=l 

where the number set {ai,...,aQ} for the asymptotic so¬ 
lution (t —>■ -boo) will depend on initial conditions (for 
the dynamical consideration). 

Let us show that the found solution m is a generaliza¬ 
tion of well-known steady-state for the constant operator 
L in Eq. ([2]): 

Lpst-st — 0 i ^tPst-st — 0 i (^j Pst-st) — 1- (1^) 

It is obviously that the stationary case can be considered 
as periodic case with an arbitrary value of the period T. 
Therefore, according to the above analysis, the state Pst-st 
also must satisfy the Eq. (ITSl) for arbitrary t, T. Indeed, 
in the case of L=const we have (for arbitrary t,T): 

^ +00 .. 

A{t + T,t) = e^^ = 1 + . (18) 

k^l 

Taking into account Eq. m, we obtain the perfect com¬ 
pliance with Eq. din]): A{t + T,t)Pst-st=Pst-st- 

Let us compare our results with the general Floquet 
theory (e.g., see [1^). On this theory the solutions of the 
differential equation system dU with periodic coefficients 
can be presented as following form: 

p{t) = D{t)e^*d, (19) 

where the matrix D{t)=D{t + T) is periodic, if is a cer¬ 
tain constant matrix connected with initial matrix L(t), 


a is an arbitrary constant vector. To have the periodic so¬ 
lution the matrix K must satisfy certain criteria. From 
this viewpoint, it is most probably that in the present 
paper we have proved the existence of the eigenvector do 
with eigenvalue 0 for the matrix K: Kdo=0. In this case, 
takin into account e^‘do=do, from Eq. ([T^ we find the 
following solution: 

p(t) = D{t)e^*ao = D(t)do , (20) 

which obviously is periodic [due to the periodicity of 
D{t)]. We believe that exactly this scenario is realized for 
density matrix equation. However, this problem requires 
an additional mathematical investigation. 

Note that some authors supposed (without proof) the 
existence of the periodic solution © for some certain 
problems. In this case they usually used the Fourier anal¬ 
ysis for numerical calculations (e.g., see [iSlil). How¬ 
ever, as far as we know, the general mathematical proof 
of the existence theorem was not previously presented in 
scientific literature. At the same time a special interest 
has the direct and simple method, which allows us to 
calculate the periodic solution on the base of Eq. (na 
without Fourier expansion, i.e. our approach automati¬ 
cally takes account of all frequency components. It can 
radically simplify the calculations for periodically-driven 
systems. 

Let us describe one of possible numerical algorithms 
based on Eq. m- We consider an arbitrary periodic de¬ 
pendence of the operator L{t) (see Fig[T|). For instance, 
under atom-field interaction such a dependence can be 
produced by the modulation of the field parameters (am¬ 
plitude, phase, polarization, etc.). The selected time in¬ 
terval [toito + T] we divide into N small subintervals, 
where tN=to + T. The character of partition (uniform 
or nonuniform discrete mesh) and number of subinter¬ 
vals are determined in conformity with studied problem. 
The dependence L{t) we will approximate by step func¬ 
tion (see Fig. dj), where the matrix L{t) has the constant 
value inside of subinterval (tm-iAm]- In this 

case the vector pitg) in initial point to is determined by 
Eq. (ITSl) . where the evolution operator A{to + T,to) has 
the form of an ordered product of the matrix exponents: 

m—N 

A{to + T,to) = E = (21) 

m—1 

The vectors p{tm) in other points of the interval [to, to+T] 
are determined by the recurrence relation: 

p{tm) = )p(t„_i). (22) 

Below we give two examples, which demonstrate the 
simplicity and high efficiency of our method. First ex¬ 
ample is the so-called dark resonance, which is formed 
in three-level A-system by bichromatic field (see Fig.d]). 
















4 


This resonance occurs when a difference between the op¬ 
tical frequencies (wi — UJ 2 ) is varied near the transition 
frequency between the lower energy levels |1) and |2): 
coi — C 02 ~ A. Currently, such resonances at the hyper- 
fine structure of the alkali metal atoms (first of all, for 
®^Rb and ^^^Cs) underlie the chip-scale atomic clocks, 
which have a great practical importance. At the same 
time the harmonic phase modulation (frequency devia¬ 
tion) of miniature semiconductor lasers (so-called VIC- 
SEL) is ordinarily used in practice. In this case the field 
can be written as 

E{t) = -h c.c.; ip{t) = A sm(iyt ), (23) 


| 3 > 



FIG. 2: Three-level system where the dark resonance can be 
formed. 


where Eq is the field amplitude, uj is the central frequency 
of laser, A and u are the amplitude and frequency of 
phase modulation, respectively. It is obvious that this 
field is polychromatic, in which the frequency difference 
between the adjacent components equals to i'. The phase 
modulation at full- or half of the hyperfine splitting is 
usually used: or A/2. However, the dark reso¬ 

nance can be observed in the more general case of vk.IS./I 

{ 1 = 1 , 2, 3 ,...). , . . . 

For theoretical description of this problem the Fourier 
expansion of the field (1^51) is traditionally used. In this 
case, to solve the density matrix equation one can ex¬ 
tract two resonant components (e.g., see (T^ ), which 
are involved in the absorption and form the dark reso¬ 
nance per se. While the remaining frequency components 
are taken into account only from the viewpoint of field- 
induced shifts of the clock transition. It is evident, that 
such approach is quite cumbersome and incomplete. In 
our method the Fourier expansion is not used and we di¬ 
rectly calculate the time-dependent periodic p{t). In the 
resonance approximation the equations for the density 
matrix components p^f. {t) have the following form: 

-I- F - iSi - i(p'{t)]p 3 i = illiipii - P33) + i^2P2i 
[ 5 t -I- F - iS2 - i(p'{t)]p 32 = i^2{P22 - P33) + i^iPl2 
[dt -I- Go — *A]pi2 = j(H*p 32 — P13O2) ( 24 ) 

[dt -\- Go]pii = 7P33/2 -|- GoTr{p }/2 -|- i{Q,\p3i — P13H3) 
[dt + Gq\p 22 = 7P33/2 -I- GoTr{p }/2 -|- i{^ 2 P 32 — P 23 ^ 2 ) 
[dt -I- Go -I- 7]P33 = *(^^lPl3 — P31^l) + i{^2P23 — P32^2) 
P,k = Pkj 0/ k = \,2, 3); Tr{p} = pn -b P22 + P33 = 1 - 

Here i 5 i, 2 =(w — ^re one-photon detunings, 

^i=d 3 iEo/h and Q. 2 =d 32 EQ/h are the Rabi frequencies 
(dsi and ^32 are the reduced matrix elements of the 
dipole moment for the respective transitions), 7 is de¬ 
cay rate of excited level |3), F is the total decoher¬ 
ence rate (spontaneous, collision, time-of-flight) of op¬ 
tical transitions | 1 )—>-13) and |2)^-|3), Gq is relaxation 
rate of lower energy levels to the equilibrium isotropic 
state: (|1)(1| + |2)(2|)/2. Note, if we delete the deriva¬ 
tive by phase modulation :p'(t), then the equations (l24l) 
formally coincide with the description of interaction be¬ 
tween monochromatic field (at the frequency w) and 
three-level A-system. 



(v-A/2)/y,p 



FIG. 3: a) The dark resonance lineshape in the case when v is 
varied near A/2 for different values of the phase modulation 
amplitude A: A=1.0 (green dashed line) and A=2.4 (solid 
black line); b) The dependance of the dark resonance position 
(relative to A/2) on parameter A. 

Galculations are done for Hi/7=H2/7=1; Go/7=1.5x10“®; 
r/7=50; A/7=1000, 52=-(5 i=A/2. 


In particular, we have calculated by our method the 
resonance lineshapes as dependencies on the frequency 
of phase modulation j/RiA/2 for different values of ampli¬ 
tude A (see Fig. [3^). The taken parameters are typical for 
atomic cells with vapor of ®^Rb in buffer gas. Fig. [5 ]d dis¬ 
plays the dependence of field shift of dark resonance (5shift 
(with respect to A/2) on modulation amplitude A. There 
is the value set of A, for which the field shift vanishes (see 
also [ 2 ^). In Fig. 2] the similar graphics are shown for 
the case of phase modulation frequency ^RiA/4. From 
Fig. 2^ it is seen that the amplitude of dark resonance 
can be comparable with the case of i^ssA/2 (compare 
with Fig. |3^). 
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a) 


(v-A/4)/y,p 



FIG. 4: a) The dark resonance lineshape in the case when v is 
varied near A/4 for the phase modulation amplitude A=3.59; 
b) The dependance of the dark resonance position (relative 
to A/4) on parameter A. 

Calculations are done for 111 / 7 = 02 / 7 = 1 ; Go/7=1.5x 10“®; 

r/7=50; A/7=1000 , (52=-5i=A/2. 


To show other examples, let us consider two-level sys¬ 
tem (with unperturbed frequency Wq) interacting with 
the field at the frequency w: i?(t)=Re{f (t) 
where the amplitude modulation £(t) has the form of 
rectangular-pulse periodic sequence (see Fig. [SJi), which 
can be considered as limiting case of the Ramsey spec¬ 
troscopy. In this case the field has an equidistant spec¬ 
trum + mfr (to= 0 ,± 1 ,± 2 ...) (so-called frequency 

comb), where /r=27r/T is the repetition frequency, and 
spectrum width (with respect to the central frequency lu) 
is determined by the pulse duration as l/r. The dynamic 
equations for density matrix components are given by: 

[(9t -I- r - iS\p2i = fO(t)(pii - /022)/2 
[i9t + F -I- iS\pi2 = —iO*(t)(pii — P22)l‘^ 

[dt + 7]P22 = i[^it)pi2 - p2ii^*(t)]/2 (25) 

dtPll = 7P22 + i[^*it)p21 — Pl2^it)]/‘2 
Trjp} = pii + P 22 = 1, 

where Q.{t)=d£{t)/h is the time-dependent Rabi fre¬ 
quency {d is the dipole matrix element), 5=(w—W q) is the 
detuning of the central frequency w. The constants 7 and 
F describe the relaxation rates of the exited state popu¬ 
lation (P 22 ) and optical coherence (pi 2 , ^ 12 ), respectively. 
Note that for pure spontaneous relaxation the condition 
r= 7/2 takes place. 

The most interest has a regime, in which the rep- 


Q(f) 


t 

b) 


t 


FIG. 5: Two-level system interacting with: a) periodic se¬ 
quence of rectangular pulses; b) periodic sequence of Ttt phase 
jumps. 
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FIG. 6: The lineshape for direct frequency comb spec¬ 

troscopy (see Fig.[S^) over large interval of detuning 5, where 
the narrow resonances are visually inseparable (dark areas). 
Calculations are done for r=7/2=0.02/r, r/r=0.001 and for 
different pulse area: a) fIo'7"=7r/2; b) flo'r=7r; c) IIoT=27r. 


etition frequency significantly exceeds the relaxation 
rate of atomic system: r, 7 <C/r-, and also the realiza¬ 
tion of the condition The Figs. [ 6 ^-c demon¬ 

strate the dependences of period-averaged value {p22)t 
as a function of detuning 5 for different pulse area 
Q.qt='k 12^71 ^2^. The structure of these lineshapes con- 
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FIG. 7: The lineshapes for direct frequency comb spec¬ 

troscopy (see Fig. Or) over two small intervals of detuning 
5, where the individual Ramsey resonances are visually sepa¬ 
rated. Calculations are done for: r=7/2=0.02/r, t/T=0.001, 
Q,qt='k/2. 

sist of large number of resonances at the frequencies 
uj=uJo+mfr (m=0,±l,±2,...). Their amplitudes lie be¬ 
tween two envelope curves. The dark areas in Fig. [S] 
correspond to the resonances, which are visually insep¬ 
arable for given scale. More detail picture using small 
frequency intervals is presented in Fig. [3 where the in¬ 
dividual resonances are visible. We have found that the 
lower envelope has a typical width of the order of 2 r 2 o 
(see in Fig. [ID and it looks similar to the lineshape on <5 
in the case of single-pulse Rabi spectroscopy. The up¬ 
per envelope is much more broadened (with width of 
the order of Ho/r/T^f^o in the case of r2oT>l). More¬ 
over, this envelope curve has a series of relatively nar¬ 
row resonance-like dips (practically to zero), which can 
be characterized as quasi-forbidden spectroscopic zones. 
The positions of these zones are determined from the con¬ 
dition: Ti/rig -I- 6'^=2k'K (fc=l,2,...), and their widths are 
proportional to 7 . 

Note the high efficiency and simplicity of our method in 
the case of periodic pulse modulation. Indeed, for above 
considered task with rectangular pulses (see Fig. [S^) the 
operator A{tQ + Tjtg) includes the product of only two 
matrix exponents: 

A{to + T, to) = g(T-r)i:(n=o) ^ 1 ^ 26 ) 

In result, the calculation of the dynamic steady-state 
in concordance with dig becomes quite simple. More¬ 
over, for two-level atom we have obtained the analytical 
expressions (exact and approximative) for the operator 
A{to + T, to) and vector p{to). In particular, for the short 
pulses (r<Cl/ 7 , t<CT) the following approximation takes 


place: 

P 22 {to) ~ fIoe'’'’^(e^'"’^ - I) sin^(nT/2)x 
{(g( 2 r-H 7 )r _ 2 e^'^(e'^'^ - l)n 6 sm{ 6 T) sin(flr)- 

— 1) cos(5T)[no + cos(nT)]-|- 

(e^^ - e^^^)[nl cosiVtT) + S^]}-\ (27) 

where 'T=(T — t) is a free evolution time between pulses, 
-I- is the generalized Rabi frequency. For the 
pure spontaneous relaxation (r= 7 / 2 ) the expression (l?7l) 
can be transformed to the form 

P22{to) ~ Oq exp(7T/2) sin^(flr/2)x 

{212^ cosh( 77 ”/2) -I- 212(5 sm{ 6 T) sin(r2T) — 
008(^7”) [fig + + ^ 0 ) cos(r 2 T)]}“^. (28) 

In addition, under condition of t<^T we have a close 
approximation: 

{p22)T«P22(to)[l-e-^^]/7T (29) 

with an absolute accuracy of small value (r/T)<Cl. Thus, 
to reproduce with a good accuracy the lineshape of {p 22 )t 
we can calculate only the dependence ^ 22 (^ 0 ) on 5. Note 
that the Eq. (12^ can be used for an arbitrary pulse form 
(not only rectangular). 

In contrast, for frequency comb spectroscopy to nu¬ 
merically solve the Eq. @ using the Fourier analy¬ 
sis we should use the following decomposition: p{t) = 
(m = 0 , ± 1 , ± 2 ,...), where the components 
Pm satisfy the certain recurrent relations. As an example, 
the calculations of the dependencies in Fig. [B] require the 
using of more than 10® frequency components (for field 
and density matrix). It leads to the huge computational 
burden (comparable with supercomputer). 

The last numerical example is shown in Fig. |8l There 
are the spectroscopy lineshapes for the field modulation, 
which is formed by ± 7 r phase jumps (see Fig. lib)- In the 
case of flo 7 , these interesting dependencies have the 
lower envelope, which has a width of the order of 2120 
(see Fig. [Hb,c,d). 

Note that for free atom gas it is necessary to take into 
account the motion of atoms, which results in dependence 
of the density matrix on velocity v. In this case the peri¬ 
odic steady-state p{t -I- T, •v)=p{t, v) should be found for 
each velocity group, and then the spectroscopic signal is 
computed with the use of velocity averaging. In addi¬ 
tion, the performed in this article consideration of tem¬ 
poral periodic modulation can be adapted for the case 
of the spatially periodic modulation of the field parame¬ 
ters (amplitude, phase, polarization, etc.). It implies the 
proof of existence and the calculation of spatially periodic 
steady-state of the atoms in the spatial lattices formed 
by electromagnetic held (including the problems of laser 
cooling and trapping of atoms in optical lattices). 

In conclusion, in the frame of density matrix formal¬ 
ism we have rigorously proven the existence theorem of 
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FIG. 8: The spectroscopy lineshapes for the field modulation, 
which is formed by ±7r phase jumps (see Fig. EId). Calcula¬ 
tions are done for T—'y/2=0.1fr and for different Rabi fre¬ 
quency: a) Qo—O-^fr', b) S7o=5/r; c) Q,o~Wfr', d) f2o=15/r. 


the periodic steady-state for an arbitrary periodically- 
driven system. Due to the relaxation processes this state 
is realized as an asymptotics {t^+oo) independently of 
initial conditions. The presented proof simultaneously 
contains a computational algorithm, which uses neither 
Floquet nor Fourier theories, i.e. our method automat¬ 
ically takes into account all frequency components. It 
radically simplifies the calculations for arbitrary types 
of periodic modulation (including the ultrashort pulses) 
and opens up great possibilities for analysis and devel¬ 
opment of new methods in the laser physics, nonlinear 
optics and spectroscopy. The importance of obtained re¬ 
sults becomes especially obvious due to the infinite va¬ 
riety of possible periodic actions for different quantum 
systems. 
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